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Abstract 

A measurement of the beauty production cross section in ep collisions at a centre-of-mass 
energy of 3 19 GeV is presented. The data were collected with the HI detector at the HERA 
collider in the years 1999-2000. Events are selected by requiring the presence of jets and 
muons in the final state. Both the long lifetime and the large mass of ^-flavoured hadrons 
are exploited to identify events containing beauty quarks. Differential cross sections are 
measured in photoproduction, with photon virtualities Q 2 < 1 GeV 2 , and in deep inelastic 
scattering, where 2 < Q 2 < 100 GeV 2 . The results are compared with perturbative QCD 
calculations to leading and next-to-leading order. The predictions are found to be somewhat 
lower than the data. 
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1 Introduction 



A measurement is presented of open beauty production ep — > ebbX in ep collisions with the 
HI detector at HERA. The measurement spans the kinematic range from the domain of photo- 
production, in which the exchanged photon is quasi-real (Q 2 ~ 0), to the region of electropro- 
duction, or deep inelastic scattering (DIS), with photon virtualities 2 < Q 2 < 100 GeV 2 . For 
beauty production, calculations in perturbative quantum chromodynamics (pQCD) are expected 
to give reliable predictions, as the mass m& of the b quark ( m& ~ 5 GeV) provides a hard scale. 
With the phase space covered in this analysis the interplay of the hard scales m&, Q 2 and the 
transverse momenta of the b quarks can be probed. 

First measurements of the beauty cross section at HERA [1,2] were higher than pQCD pre- 
dictions calculated at next-to-leading order (NLO). Similar observations were made in hadron- 
hadron collisions [3] and also in two-photon interactions [4]. Recent beauty production mea- 
surements from the HI [5] and ZEUS Collaborations [6] are in better agreement with QCD 
predictions or again somewhat higher [7]. 

In this paper, photoproduction events with at least two jets (jj) and a muon (fx) in the final 
state are used to measure the beauty cross section for ep — > ebbX — > ejjfiX'. In deep inelastic 
scattering, the process ep — > ebbX — > ejfiX' is measured with at least one jet and a muon 
in the final state. For the first time at HERA, two distinct features of 5-hadrons are exploited 
simultaneously to discriminate events containing beauty from those with only charm or light 
quarks: the large mass and the long lifetime. The £?-hadron mass leads to a broad distribution 
of the transverse momentum p\ d of decay muons relative to the beauty quark jet direction. The 
5-hadron lifetime is reflected in the large impact parameters 5 ~ 200 /im of decay muon tracks 
relative to the primary vertex. The precision measurement of muon track impact parameters is 
made possible by the HI Central Silicon Track detector [8]. The fractions of b quark events in 
the data samples are determined from a fit to the two-dimensional distribution of p\ el and 5. 

This paper is organized as follows. In section 2, an introduction to the physics of beauty 
production in ep collisions is given. The relevant features of the HI detector are described in 
section 3. Section 4 describes the event selection. The Monte Carlo simulations and NLO QCD 
calculations are presented in sections 5 and 6. Comparisons of the data samples with the Monte 
Carlo simulations are shown in section 7. The fit procedure used to determine the b- fraction and 
the systematic errors of the measurement are explained in sections 8 and 9. Finally, the results 
are presented in section 10 and conclusions are drawn in section 11. 

2 Heavy Quark Production in ep Collisions 

In pQCD, at leading order, two distinct classes of processes contribute to the production of 
beauty quarks in ep collisions at HERA. In direct-photon processes (figure la), the photon 
emitted from the positron enters the hard process 75 — > bb directly. In resolved-photon pro- 
cesses (figures lb to d), the photon fluctuates into a hadronic state before the hard interaction 
and acts as a source of partons, one of which takes part in the hard interaction. Resolved photon 
processes are expected to contribute significantly in the photoproduction region, in which the 
photon is quasi-real, and to be suppressed towards higher Q 2 . 
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Direct-7 Resolved-7 




a) 75f-Fusion b) Hadron-like c) fe-Excitation d) 6-Excitation 



Figure 1: Beauty production processes in leading order pQCD. 

For the kinematic range covered in this analysis, the majority of events are in the region 
Pt,b ^ m-b, where p t ,b is the momentum of the outgoing b quark transverse to the photon-proton 
axis in the photon-proton centre-of-mass frame. In this region, NLO calculations in the massive 
scheme [9, 10] are expected to give reliable results. In this scheme, u, d and s are the only 
active flavours in the proton and the photon, and charm and beauty are produced dynamically 
in the hard scattering. At large transverse momenta p tjb 3> m b or large Q 2 3> (2m&) 2 , the 
massive scheme becomes unreliable due to large terms in the perturbation series of the form 
a s \n(p 2 b /ml) or a s ln(Q 2 /m 2 ). In this kinematic range, the massless scheme [11] can be used, 
in which charm and beauty are treated as active flavours in both the proton and the photon, in 
addition to u, d and s. In this scheme, so-called excitation processes occur in which the beauty 
quark is a constituent of the resolved photon (sketched in figures lc and d) or of the proton. 

In this analysis the measurements are compared with NLO calculations in the massive 
scheme for both photoproduction [9] and DIS [10]. NLO calculations in the massless scheme 
are not yet available for the exclusive final state considered in this measurement. The data 
are also compared with the predictions of the Monte Carlo simulations PYTHIA [12], RAP- 
GAP [13] and CASCADE [14]. In the Monte Carlo simulations leading order matrix elements 
are implemented and higher orders are approximated using parton showers radiated from the 
initial and final state partons. PYTHIA and RAPGAP use the DGLAP [15] parton evolution 
equations, while CASCADE contains an implementation of the CCFM [16] evolution equa- 
tion. In CASCADE the direct process — > bb is implemented using off-shell matrix elements 
convoluted with A; r unintegrated parton distributions in the proton. 



3 Detector Description 

The HI detector is described in detail in [17] and only the components most relevant for this 
analysis are briefly discussed here. A right handed coordinate system is employed that has its 
z-axis pointing in the proton beam, or forward, direction and x (y) pointing in the horizontal 
(vertical) direction. Charged particles are measured in the Central Tracking Detector (CTD) 
which covers the range in pseudo-rapidity between —1.74 < rj < 1.74 1 . The CTD comprises 

'The pseudo-rapidity i] corresponding to a polar angle 9 is given by r/ = — In tan(0/2). 
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two large cylindrical Central Jet Chambers (CJCs) and two ^-chambers arranged concentrically 
around the beam-line within a solenoidal magnetic field of 1.15 T. The CTD also provides 
trigger information which is based on track segments in the r-(f> plane measured in the CJCs and 
the z -position of the interaction vertex obtained from a double layer of multi-wire proportional 
chambers. 

The CTD tracks are linked with hits in the Central Silicon Track detector (CST) [8], which 
consists of two 36 cm long concentric cylindrical layers of silicon strip detectors, surrounding 
the beam pipe at radii of 57.5 mm and 97 mm from the beam axis. The CST covers a pseudo- 
rapidity range of —1.3 < 77 < 1.3 for tracks passing through both layers. The double- sided 
silicon detectors provide resolutions of 12 fxm in r-(f> and 25 fxm in z. Average hit efficiencies 
are 97% (92%) in r-<f> (z). For a CTD track with CST r-<f> hits in both layers, the transverse 
distance of closest approach dca of the track to the nominal vertex in x-y can be measured 
with a resolution of a dca m 33 pim © 90 /xm/p t [GeV], where the first term represents the 
intrinsic resolution (including alignment uncertainties) and the second term is the contribution 
from multiple scattering in the beam pipe and the CST; p t is the transverse momentum of the 
track. 

The track detectors are surrounded in the forward and central directions, — 1.5 < rj < 3.4, by 
a fine grained Liquid Argon calorimeter (LAr) and in the backward region, —4.0 < rj < —1.4, 
by a lead- scintillating fibre calorimeter SpaCal [18] with electromagnetic and hadronic sections. 
The SpaCal is used primarily in this analysis to detect the scattered electron in DIS events and 
to select photoproduction events, in which case the scattered electron is not detected. The 
calorimeters are surrounded by the solenoidal magnet and the iron return yoke which is instru- 
mented with 16 layers of limited streamer tubes in the range —2.5 < rj < 3.4. In the central 
pseudo-rapidity range, studied in this paper, the instrumented iron allows high efficiency detec- 
tion of muon tracks with pf > 2 GeV. 

The ep luminosity is determined by measuring the QED bremsstrahlung (ep — > ep^y) event 
rate by detecting the radiated photon in a calorimeter located at z = — 103 m. 

4 Event Selection and Reconstruction 

The data were recorded in the years 1999 and 2000 and correspond to an integrated luminosity 
of ~ 50 pb -1 . During this time HERA was operated with positrons of 27.6 GeV and protons 
of 920 GeV energy. The events were triggered by requiring that there be signals from the 
central drift chambers and the multi-wire proportional chambers in coincidence with signals 
from the scattered positron in the backward calorimeter (DIS sample) or with signals from the 
instrumented iron (photoproduction). 

At least one muon is required with a transverse momentum pf > 2.5 GeV. Muons are 
identified as track segments in the barrel part of the instrumented iron. The iron track segments 
must be well matched to a track reconstructed in the CTD. At least two CST r-0-hits have to 
be associated with the muon track and it is required that the combined CTD-CST r-0-track 
fit probability exceeds 10%. The muon momentum is reconstructed using the CTD-CST track 
information. The CST hit requirements for the muon track restrict the allowed range of ep 
interactions along the z-axis to \z vtx \ < 20 cm. 
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Jets are reconstructed using the inclusive k t algorithm [19, 20] in the p t recombination 
scheme (see [21]), giving massless jets. The algorithm, with a distance parameter in the 77- 
<p plane of 1, is applied to all hadronic final state particles, which are reconstructed using a 
combination of tracks and calorimeter energy deposits [22]. The muon, as measured in the 
CTD and CST, is one of the particles which is input to the jet algorithm. The muon track is 
required to be associated with one of the selected jets by the jet algorithm. 

Photoproduction events are selected by demanding that there be no electromagnetic cluster 
in the backward and central calorimeter with an energy of more than 8 GeV. This cut restricts 
the accepted range of negative four-momentum transfer squared to Q 2 < 1 GeV 2 with a mean 
of about 0.07 GeV 2 . For the photoproduction sample, a cut on the inelasticity 0.2 < y < 0.8 is 
applied, where y is reconstructed using the relation y = J2 h (E — p z )/2E e [23]. Here, E and 
p z are the energies and ^-components of the momenta of the hadronic final state particles, h, 
and E e is the positron beam energy. The final photoproduction event sample consists of 1745 
events. The number of events containing more than one muon candidate is less than 1%. The 
jet algorithm is applied in the laboratory frame and at least two jets are required with transverse 
momenta p\ e 1(2) > 7(6) GeV in the pseudo-rapidity range \i] jet \ < 2.5. The fraction of the 
photon energy entering the hard interaction is estimated using the observable 

where the sums in the numerator run over the particles associated with the two jets and that in the 
denominator over all detected hadronic final state particles. For the direct process (figure la), 
x° bs approaches unity, as the hadronic final state consists of only the two hard jets and the proton 
remnant in the forward region which contributes little to ^ h (E — p z ). In resolved processes 
x obs can be small. 

DIS events are selected by requiring a scattered positron signal with an energy of at least 
8 GeV in the SpaCal. To suppress photoproduction background and to reduce the fraction of 
events with significant initial state QED radiation, events are rejected if ^2(E — p z ) < 45 GeV. 
Here, E — p z is summed over all final state particles including the scattered positron. The 
kinematic variables Q 2 and y are reconstructed using the eS method [24], which combines 
the hadronic final state and the positron measurements. The scaling variable x is subsequently 
calculated using the relation x = Q 2 /ys, where s is the ep centre-of-mass energy squared. 
Events are selected in the range 2 < Q 2 < 100 GeV 2 and 0.1 < y < 0.7. The jet algorithm is 
applied in the Breit frame [25] and at least one jet with transverse momentum pfj^ 1 > 6 GeV 
is required, with which the muon must be associated. The final DIS event sample consists of 
776 events. 

Table 1 summarises the selection cuts for the two samples. The selection cuts for the pho- 
toproduction sample are somewhat tighter, due to the muon trigger acceptance and to suppress 
background events. 

5 Monte Carlo Simulations 

The Monte Carlo generators PYTHIA [12] and RAPGAP [13] are used for the description of 
the signal and background distributions, the determination of efficiencies and acceptances and 
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Photoproduction 


DIS 


Q 2 [GeV 2 ] 


< 1 


2... 100 


y 


0.2 ... 0.8 


0.1 ... 0.7 


Frame 


laboratory 


Breit 


#jets 


> 2 


> 1 


[GeV] 


>7(6) 


> 6 


\Vlab\ 


< 2.5 


< 2.5 


fj, Iron link probability 


> 10% 


>5% 


fi CST link probability 


> 10% 


> 10% 


# CST hits 


> 2 


> 2 




-0.55 ... 1.1 


-0.75 ... 1.15 


Pt [GeV] 


> 2.5 


> 2.5 


# events 


1745 


776 



Table 1 : Selection cuts for the photoproduction and DIS data samples and the number of se- 
lected events. 



for systematic studies. The track resolutions were adjusted to describe the data. In addition, sys- 
tematic cross checks and estimates of model dependence are performed using the CASCADE 
generator [14]. The measured beauty production cross sections are also compared with the pre- 
dictions of these three generators. The basic parameter choices for the various pQCD programs 
are summarised in table 2. 

All three Monte Carlo generators are used to produce large samples of beauty and charm 
events with decays into muons ep — > ebbX —> fiX' or ep eccX — > fiX'. The Peterson 
fragmentation function [26] is used for the hadronisation of the heavy quarks. For systematic 
cross checks, samples using the Lund string fragmentation model [27] are generated. Each of 
these Monte Carlo samples corresponds to at least forty times the luminosity of the data. In 
addition, PYTHIA and RAPGAP event samples for light quark (uds), c and b events without 
muon requirements are generated with six times the luminosity of the data. The latter samples 
are used for the simulation of the background due to hadrons misidentified as muons and decays 
of light hadrons into muons. All generated events are passed through a detailed simulation of 
the detector response based on the GEANT program [28] and reconstructed using the same 
reconstruction software as used for the data. 

For the measurements in photoproduction, PYTHIA is used in an inclusive mode in which 
direct and resolved events are generated using massless matrix elements for all quark flavours 
(MSTP(14)=30 [12]). Beauty (charm) events are separated from light quark events by requiring 
that there be at least one beauty (charm) quark in the list of outgoing hard partons. Approx- 
imately 35% of the PYTHIA beauty cross section in the measured range is due to resolved 
photon processes and these are dominated by the flavour excitation component. For the mea- 
surements in the DIS region, RAPGAP is used to generate the direct production process in the 
massive mode (IPRO=14 [13]). RAPGAP is interfaced with the program HERACLES [29] 
which simulates QED initial and final state radiation. Additional event samples are generated 
using the Monte Carlo generator CASCADE [14]. Comparisons of the measurements with the 
predictions of CASCADE are made using version 1.2 of the program and the J2003 parton 
density distributions [30-32]. 
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PYTHIA 


RAPGAP 


CASCADE 


FMNR 


HVQDIS 


Version 


6.1 


2.8 


1.00/09; 1.2 




1.4 


Proton PDF 


CTEQ5L [33] 


CTEQ5L 


JS2001 [14] 


CTEQ5M [33] 


CTEQ5F4 [33] 








J2003 [30-32] 






Photon PDF 


GRV-G LO [34] 






GRV-G HO [34] 




^q ] cd [GeV] 


0.192 


0.192 


0.2 


0.326 


0.309 


Renorm. scale [/% 


m o + Ptaa 


Q 2 +Pm 




m l + PthT, 

o 1 too 


m l + Pthb 

o 1 tOO 


Factor, scale u 2 ( 


m?, 4- y>+„x 

q 1 rtqq 


Q 2 + Ptnn 


s + Qf 


ml + V^r 
o ^tbb 


ml + P^r 
b 1 tbb 


m b [GeV] 


4.75 


4.75 


4.75 


4.75 


4.75 


m c [GeV] 


1.5 


1.5 


1.5 






Peterson e& 


0.0069 


0.0069 


0.0069 


0.0033 


0.0033 


Peterson e c 


0.058 


0.058 


0.058 







Table 2: Parameters used in the leading order Monte Carlo simulations and the NLO programs. 
Here ji r and fif denote the renormalisation and factorisation scales, m q the heavy quark masses, 
Pt q q the average of the transverse momenta of the two heavy quarks, s and Q\ the heavy quark 
system centre-of-mass energy squared and transverse momentum squared, respectively, and e q 
the Peterson fragmentation parameters. 



All generators use the JETSET part of the PYTHIA [12] program to simulate the hadronisa- 
tion and decay processes. The branching ratios for the direct semileptonic decays b — > jiX and 
for the indirect decays into muons via charm, anticharm, r and J/ ^ decays are in agreement 
with the world average values [35]. The total branching ratio for beauty decays into muons is 
21% [35]. The decay lifetimes of the beauty and charm hadrons are set to the values reported 
in [35]. The muon momentum spectrum in the rest frame of the decaying b- flavoured hadrons, 
as modeled by JETSET, is in agreement with the spectrum measured at e + e~ colliders [36, 37]. 

6 NLO QCD Calculations 

The NLO pQCD calculations are performed in the massive scheme using the program FMNR [9] 
in the photoproduction regime and the program HVQDIS [10] for the DIS case. Both programs 
provide weighted parton level events with two or three outgoing partons, i.e. a b quark, a b quark 
and possibly an additional light parton. The calculations are performed in the MS-scheme using 
the parameters given in table 2. 

The b quark is 'hadronised' into a ^-flavoured hadron by rescaling the three-momentum of 
the quark using the Peterson fragmentation function [26] with the parameter e b = 0.0033 [38]. 
The programs are extended to include the decay of the b- flavoured hadron into a final state with 
a muon. The muon decay spectrum is taken from JETSET [12] and includes direct and indirect 
decays of ^-flavoured hadrons into muons. Parton level jets are reconstructed by applying the 
k t jet algorithm to the outgoing partons. 

Corrections to the hadron level are calculated using the PYTHIA and RAPGAP Monte Carlo 
event generators. PYTHIA and RAPGAP parton level jets are reconstructed from the generated 
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quarks and gluons after the parton showering step. At the hadron level, jets are reconstructed 
by applying the jet algorithm to all final state particles, after the decay of the beauty or charm 
hadrons. The jet and the muon selection cuts are applied to the generator samples. In each kine- 
matic bin of the measurement, the ratio of the PYTHIA or RAPGAP hadron level and parton 
level cross sections is calculated and applied as a correction factor to the NLO calculation. The 
parton to hadron level corrections range typically from —30% to +5% in both photoproduction 
and DIS. The corrections are negative at small muon and jet transverse momenta and positive at 
the largest transverse momenta. The corrections obtained using CASCADE are consistent with 
the values from PYTHIA (photoproduction) and RAPGAP (DIS). 

The theoretical uncertainties of the NLO calculations are estimated in the following way: 
The b quark mass and the renormalisation and factorisation scales are varied simultaneously 
from m b = 4.5 GeV and ji r — /// = m T /2 to m b = 5 GeV and fi r — /if — 2m T , where 



m T = y ml + p 2 bb - , and p tbb is the average of the transverse momenta of the two b quarks. This 
leads to a maximum change of the cross section of typically 25% in photoproduction (FMNR) 
and 15-20% in DIS (HVQDIS). The cross section variation when using other proton structure 
functions such as MRSG or MRST1 [39] is less than 8% in all regions of the measurement. 
The uncertainty due to variations of the fragmentation parameter e b by 25% is below 3%. These 
cross section variations are added in quadrature to estimate the total systematic uncertainty of 
the NLO predictions for each bin of the measurement. 



7 Comparison of the Data with Monte Carlo Simulations 



Detailed comparisons are performed of the data with the Monte Carlo simulations. The good 
agreement of the PYTHIA simulation with the photoproduction data is illustrated in figure 2. 
The distributions of the muon transverse momentum p%, the pseudo-rapidity rf, the jet trans- 
verse momentum p t e 1(2) and x° bs are shown. The data are compared with the sum of the con- 
tributions from beauty, charm and light quark events, the relative fractions of which are taken 
from the two-dimensional fit discussed in section 8. The number of events in the simulation 
is normalised to that of the data. It is observed that the shapes of the distributions are rather 
similar for the beauty, charm and light quark events except for the muon transverse momentum, 
where the spectrum is harder for beauty than for the other components. In figure 3, the distri- 
butions for the DIS sample are shown. The photon virtuality Q 2 , the inelasticity y, Bjorken-rc, 
the muon transverse momentum p% and the transverse momentum pf^lf of the selected jet in 
the Breit frame are well described by the RAPGAP Monte Carlo simulation. The CASCADE 
Monte Carlo simulation also provides a good description of the data in both photoproduction 
and DIS (not shown). 



8 Determination of Beauty and Background Contributions 



The signed impact parameter 5 of the muon track and the transverse momentum p r t el of the 
muon track relative to the axis of the associated jet are used to determine the fraction of beauty 
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events in the data. For each muon candidate, 5 is calculated in the plane transverse to the beam 
axis. The magnitude of 5 is given by the dca of the track to the primary vertex. The sign is 
defined as positive if the angle between the jet axis and the line joining the primary vertex to 
the point of closest approach of the track is less than 90°, and is defined as negative otherwise. 
Figures 4a and 5a show the distributions of 5 for the photoproduction and DIS samples, respec- 
tively. The decay of long-lived particles mainly leads to positive impact parameters, whereas 
particles produced at the vertex yield a symmetric distribution centered at zero with finite width 
due to the track and primary vertex resolutions. 

The transverse beam interaction region at HERA, in the following termed the 'beam spot' , 
has an approximately Gaussian profile. For the data period studied in this paper the beam spot 
size is determined to be about 145 /im in the horizontal and 25 fim in the vertical direction. 
For each event the ep collision point is determined in a primary vertex fit using the weighted 
average of the beam spot position, taking the above widths as errors, and position information 
from selected tracks in the event. The muon track under consideration is excluded from the 
fit. An average muon impact parameter resolution of 80 /im is achieved, with comparable 
contributions from the muon track resolution and the uncertainty in the primary vertex position. 

The transverse momentum p r t el of the muon track (figures 4b and 5b) is calculated relative 
to the direction of the associated jet according to the formula 

rel _ fa X (Pjet-jV)l 
Pt ~ I -> —* | 

\Pjet — Pfi\ 

The quantities and pj et are the momentum vectors of the muon and the jet in the laboratory 
frame, respectively. 

The fraction of muons in the data that originate from beauty events is determined using a 
likelihood fit to the two-dimensional distribution of 5 and p r t el in the range —0.1 < 5 < 0.1 cm 
and < p\ el < 3.6 GeV. The combination of the two independent observables 5 and p r t el in 
the fit results in a significant improvement in the statistical precision of the measurement and 
a reduced sensitivity of the measured b fraction to systematic uncertainties in the modelling of 
charm and light quark contributions. The shapes of the distributions 5 and p r t el for beauty 2 , 
charm and light quark events are obtained from the Monte Carlo simulation. In the fit, the 
relative fractions of the three components are adjusted such that the likelihood is maximised. 
The normalisation of the sum of the three components is fixed to match the data. The results of 
the fits are illustrated in figures 4 and 5. The contributions from beauty, charm and light quark 
events, with respective fractions of typically 30%, 50% and 20%, are indicated. In all bins of 
the measurement, the data are well described by the sum of the three contributions. At large 
positive values of 5 and at large values of p r t el , the beauty component (dashed line) becomes 
dominant. 

Variations of the fit procedure are investigated and are found to give consistent results. For 
example, the S and p r t el distributions are fitted independently of each other. In another fit the 
relative contributions of charm and light quark events are fixed to the predictions of the Monte 
Carlo simulation. Furthermore, the fits to the two-dimensional data distributions of 5 and p r t el 
are investigated in a beauty enriched subset of the 5 and p r t el phase space, as illustrated in the 

2 Both direct decays b — > fiX and indirect decays, e.g. b — ► cX' — ► fj,X, are taken into account. 
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figures 6 and 7. Here, the distributions of the impact parameter 5 (figures 6a and 7a) and the rel- 
ative transverse muon momentum p\ el (figures 6b and 7b) are shown for the cuts p r t el > 1.2 GeV 
and 5 > 0.01 cm, respectively. The lines show the predictions for the different contributions in 
the restricted samples when using the results of the fits to the complete samples. The expected 
enhancement of the beauty contribution is observed and the quality of the description of the 
data illustrates the consistency between the results obtained using the two observables indepen- 
dently. 

In each kinematic bin of the measurement, the fit to the two-dimensional data distribution of 
5 and p[ el is performed, using the data and Monte Carlo samples in that bin. The fitted number 
of muons coming from beauty events is translated into a cross section by correcting for detector 
efficiencies, acceptances and radiative effects and by dividing by the integrated luminosity. 



9 Systematic Uncertainties 

Systematic uncertainties of the cross section measurement are evaluated by variations applied 
to the Monte Carlo simulations. The dominant errors come from the muon identification and 
muon track linking efficiencies, the modelling of the resolution of the muon impact parameter 
and the fragmentation models. The systematic errors assigned to the measured cross sections 
are listed in table 3. 

The muon track reconstruction efficiency in the CTD is known to a precision of about 2%. 
An additional uncertainty of 2% comes from the requirement that two CST hits be associated 
with the central track, yielding a total uncertainty for the track reconstruction efficiency of 
3%. The uncertainty of the muon identification efficiency, including the reconstruction in the 
instrumented iron and the linking with the central track, is about 5%. 

The systematic error arising from the uncertainties of the CTD and CST track resolutions is 
estimated by varying the muon impact parameter resolution in the Monte Carlo simulations by 
10%. This leads to cross section changes of 7%. To substantiate this result the following cross 
checks are performed and found to be consistent. The core and the tails of the distribution of 
the impact parameter resolution are varied separately. The description of the beam spot ellipse 
is tested by determining the cross sections separately for two independent samples with either 
more horizontal or more vertical muons. In addition, the muon impact parameter is calculated 
with respect to the average ep collision point instead of the primary vertex. 

The reconstruction of the direction of the jet associated to the muon is studied by varying 
the resolutions of the reconstructed jet directions in the Monte-Carlo simulations. The effect 
on the measured cross sections is about 2%. The jet energy scale uncertainties are estimated by 
varying the LAr energy scale in the Monte Carlo simulations by 4%. This leads to cross section 
changes of up to 4%. 

The trigger efficiencies are 70 ± 3% for the photoproduction sample and 85 ± 3% for the 
DIS sample, respectively. For the DIS sample, the uncertainty associated with reconstruction 
and identification of the scattered positron is estimated to be less than 2%. The luminosity 
measurement contributes a global 1.5% error. 
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The dependence on the physics model used for the beauty signal and the charm background 
is studied using the CASCADE Monte Carlo generator instead of PYTHIA or RAPGAP, leading 
to cross section variations of about 5%. Using the Lund [27] fragmentation model instead of 
the Peterson fragmentation function [26] causes changes in the measured cross sections of up 
to 7%. 

The modelling of the decays of the 6-flavoured hadrons has been tested by varying the 
lifetimes and branching ratios of the different hadrons within the uncertainties of the world 
average values. The effects on the measured cross section are at the 2% level. Muons from 
7r ± or K ± decays within the beam pipe and inside the sensitive volume of the CTD and CST 
exhibit a broad 5 distribution. The contribution from these events in the light quark Monte Carlo 
simulation is varied by a factor of two, leading to cross section changes of 2%. 

The above systematic studies are performed separately for each bin of the cross section 
measurement. The systematic errors are found to be of similar size for all bins. For each bin 
a total systematic uncertainty of 14% is estimated by adding all contributions to the systematic 
error in quadrature. 



Source 


Photoproduction 

Aa/a [%} 


DIS 

Aa/a [%} 


Detector efficiencies 






- Scattered positron 




2 


- Trigger efficiency 


4 


3 


- Muon identification 


5 


5 


- CST+CTD tracks 


3 


3 


- Luminosity 


1.5 


1.5 


Track reconstruction 






- 5 resolution 


7 


7 


Jet reconstruction 






- Jet axis 


2 


2 


- Hadronic energy scale 


4 


4 


MC model uncertainties: 






- PYTHIA vs. CASCADE 


5 




- RAPGAP vs. CASCADE 




5 


- Fragmentation (Peterson vs. Lund) 


7 


7 


- Fragm. fractions, BRs, lifetimes 


2 


2 


- K,n decays 


2 


2 


Total 


14 


14 



Table 3: List of systematic uncertainties as discussed in section 9. The total systematic error is 
obtained by adding all contributions in quadrature. 
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10 Results 



Differential beauty production cross sections are determined separately for the photoproduction 
and electroproduction samples. The results are listed in tables 4, 5 and 6 and displayed in 
figures 8 to 11. 

10.1 Photoproduction Measurement 

The visible range for the measurement for beauty photoproduction in dijet muon events is 
Q 2 < 1 GeV 2 , 0.2 < y < 0.8, pj* > 2.5 GeV, -0.55 < rf < 1.1, pf * 1(2) > 7(6) GeV and 
\r) jet \ < 2.5. For this range the total cross section is measured to be 

(?vis(ep -> ebbX -> ejjfiX') = 38.4 ± 3A(stat.) ± 5A(sys.) pb. 

The NLO QCD calculation performed in the massive scheme with the FMNR program [9], as 
described in section 6, yields for the same kinematic range a value of 23.8lg ;^pb which is 1.5 
standard deviations below the data. The Monte Carlo programs PYTHIA and CASCADE also 
predict a lower cross section than that measured in the data (see table 6). The results of all three 
calculations are in good agreement with each other. An analysis using an independent HI data 
sample was performed in [40] giving results consistent with this measurement. 

Differential cross sections for beauty production are measured as a function of several kine- 
matic variables, shown in figure 8 and listed in table 4. The bins in which the measurement is 
made are identical to the bins in which the theory curves are presented. The measured cross 
sections are quoted at the point in the bin at which the bin-averaged cross section equals the 
differential cross section, according to the Monte Carlo simulation. The data are compared 
with the expectations of the FMNR NLO QCD calculation and the PYTHIA and CASCADE 
generators. 

The differential cross section measured as a function of the muon pseudo-rapidity r/^ (fig- 
ure 8a) is flat in the phase space covered. The NLO QCD calculation describes the shape well. 
This is also true for both PYTHIA and CASCADE. The measurement agrees well with the val- 
ues obtained by the ZEUS experiment [6] in their two central muon pseudo-rapidity bins, which 
cover a similar phase space. 

The differential cross sections measured as a function of the muon transverse momentum 
Pt and of the transverse momentum of the leading jet jP t et (figure 8b and c) fall steeply with 
increasing transverse momentum. The NLO calculation clearly predicts a less steep behaviour 
and is lower than the data in the lowest momentum bin by roughly a factor of 2.5. At higher 
transverse momenta better agreement is observed. Similar conclusions can be drawn for both the 
PYTHIA and CASCADE predictions, although the latter predicts a slightly harder p\ et spectrum 
than the other calculations. 

Figure 8d shows the differential cross sections as a function of x° bs . A significant fraction of 
the data is found at x° bs < 0.75, i.e. in the region in which resolved photon processes (figures lb 
to d) are enhanced. In this observable the NLO calculation suffers from large uncertainties due 
to the scale variations. Furthermore, the parton to hadron level corrections are large due to the 
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fact that a single parton can produce more than one jet at the hadron level leading to migrations 
in x° bs . Within the large uncertainties, the NLO calculation describes the x° bs differential cross 
sections reasonably well. The PYTHIA simulation includes a 35% contribution from resolved 
photon processes, which are dominated by flavour excitation processes such as those shown in 
figures lc and d. Due to the large fraction of resolved photon processes, PYTHIA predicts a 
relatively high cross section value in the lowest x obs bin (figure 8d), which matches the data quite 
well. However, in the largest bin, x obs > 0.75, the PYTHIA prediction is too low. In contrast, 
CASCADE succeeds in describing the cross section of the data at large values of x° bs while it 
is too low at smaller values of x° bs . 

The results of this analysis are compared with the previous HI measurement [1] in which 
somewhat softer jet and muon cuts than in this analysis were used for the event selection. The 
measured cross section for the process ep — > ebbX — > ejjfxX' is extrapolated to the inclusive b 
quark cross section, ep — > bbX — > /iX', in the kinematic region Q 2 < 1 GeV 2 , 0.1 < y < 0.8, 
Pt > 2 GeV and 35° < 9^ < 130°, as in [1]. The extrapolation is performed using the Monte 
Carlo program AROMA [41] which was also used in [1]. The result, scaled to 820 GeV proton 
beam energy, is 107.3 ± 9.5(stat.) ± 15.1(sys.) pb, which is 2.3 standard deviations lower than 
the value of 176 ± lQ(stat.) ±™(sys.) pb obtained in [1]. 



10.2 Electroproduction Measurement 

The beauty electroproduction cross section is measured in the visible range 2 < Q 2 < 100 GeV 2 , 
0.1 < y < 0.7, pt > 2.5 GeV, -0.75 < rf < 1.15, pg™* > 6 GeV and \^ et \ < 2.5, yielding 

&vis(ep -> ebbX -> ejfiX') = 16.3 ± 2.0(stat.) ± 2.3(sys.) pb. 

The prediction of the NLO QCD calculation in the massive scheme using the program HVQDIS 
is 9.0^ 2 'g pb, which is 1.8 standard deviations below the data. The Monte Carlo programs 
RAPGAP and CASCADE also predict a lower cross section than that measured in the data (see 
table 6). 

Differential cross section measurements are presented in figures 9 and 10 and in table 5. 
The data are compared with the expectations of the HVQDIS NLO QCD calculation and the 
RAPGAP and CASCADE generators. The differential cross section as a function of the photon 
virtuality Q 2 is shown in figure 9a. The NLO calculation describes the shape well, but lies 
below the data. The prediction of CASCADE is similar to that of the NLO calculation while 
RAPGAP, which only contains the direct photon contribution to the cross section, is somewhat 
further below the data. The differential cross section as a function of the scaling variable x is 
shown in figure 9b. The various calculations also describe the shape of the data well, while the 
overall normalisation is again too low. 

In figure 10, the differential cross sections are presented as functions of the muon and lead- 
ing jet kinematics. The differential cross section measured in bins of the muon pseudo-rapidity 
r]^ (figure 10a) exhibits a rise towards the forward region, which is not reproduced by the NLO 
and Monte Carlo calculations. The differential cross sections measured as a function of the 
transverse momenta of the muon and of the jet in the Breit frame pfj^f (figures 10b and c) 
show a steep distribution, as is the case in photoproduction (figures 8b and c). The shapes of 
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the NLO QCD, RAPGAP and CASCADE predictions are all very similar. As in photoproduc- 
tion, the measured electroproduction cross sections as a function of the muon and jet transverse 
momenta show a steeper behaviour than the predictions of the NLO calculations and the Monte 
Carlo simulations and significantly exceed the predictions in the lowest bins. 

Figure 11 presents a summary of recent HERA beauty cross section measurements as a 
function of the photon virtuality Q 2 . The figure shows the ratios of the measured cross sec- 
tions [5-7] and the corresponding next-to-leading order predictions where FMNR is used for 
the photoproduction and HVQDIS for the DIS region. The dotted lines indicate the typical the- 
ory error due to scale uncertainties. General agreement is seen between the results from HI and 
ZEUS, the data tending to be somewhat above the NLO predictions. 

11 Conclusions 

Differential beauty production cross sections are measured in ep collisions at HERA both in 
photoproduction (Q 2 < 1 GeV 2 ) and in electroproduction (2 < Q 2 < 100 GeV 2 ). The event 
selection requires the presence of at least one jet (two jets) in the DIS (photoproduction) sample 
and a muon in the central pseudo-rapidity range. For the first time at HERA, beauty events are 
identified using both the transverse momentum of the muon relative to the jet axis and the large 
impact parameter of the muon. The cross sections presented here are in general agreement with 
those obtained by the ZEUS experiment. The data are compared with predictions based on NLO 
QCD calculations in the massive scheme and with the expectations of Monte Carlo generators 
which use leading order matrix elements and parton showers. The predictions from all these 
calculations are similar in both normalisation and shape. 

In both photoproduction and DIS, the total cross section measurements are somewhat higher 
than the predictions. The excess is observed mainly at small muon and jet transverse momenta, 
while at larger momenta a reasonable description is obtained. In photoproduction a significant 
contribution to the cross section is observed in the region of small values of the observable 
x° bs , where contributions from resolved photon events are enhanced. In this region the best de- 
scription of the data is given by the PYTHIA simulation, which incorporates flavour excitation 
processes in which the beauty quark is a constituent of the resolved photon or the proton. In 
DIS, the observed excess is pronounced at large muon pseudo-rapidities. The shape of the cross 
section as function of the photon virtuality Q 2 is reproduced by the QCD calculations over the 
full range covered by the measurement presented in this paper, from quasi-real photons up to 
virtualities of about 4m?. 
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Measurement 


Experimental Errors 


^-range 


tf" 


dajdrf 


stat. 


syst. 


total 






[pb] 


[pb] 


[pb] 


[pb] 


-0.55 -0.15 


-0.35 


19.1 


3.4 


2.7 


4.3 


-0.15 0.25 


0.05 


01 A 

15 .4 


4.0 


3.3 


5.2 


0.25 0.65 


0.45 


01 Q 


3.9 


3.3 


5.1 


0.65 1.10 


0.85 


O 1 o 

21.8 


3.8 


3.0 


4.9 


Pt -range 


Pt 


da/dpt 


stat. 


syst. 


total 


[GeV] 


[GeV] 


[pb/GeV] 


[pb/GeV] 


[pb/GeV] 


[pb/GeV] 


2.5 3.3 


2.9 


24.4 


3.3 


3.4 


4.8 


3.3 5.0 


4.1 


8.1 


1.1 


1.1 


1.6 


5.0 12.0 


7.2 


1.15 


0.18 


0.16 


0.24 


p\ et -range 


jet 
Pt 


dcr/dpl et 


stat. 


syst. 


total 


[GeV] 


[GeV] 


[pb/GeV] 


[pb/GeV] 


[pb/GeV] 


[pb/GeV] 


7.0 10.0 


8.8 


6.3 


0.8 


0.9 


1.2 


10.0 14.0 


11.7 


2.9 


0.4 


0.4 


0.6 


14.0 25.0 


18.3 


0.83 


0.14 


0.12 


0.18 


x° bs -range 




da/dx° bs 


stat. 


syst. 


total 






[pb] 


[pb] 


[pb] 


[pb] 


0.20 0.50 


0.35 


17.2 


4.5 


2.4 


5.1 


0.50 0.75 


0.63 


21.4 


5.2 


3.0 


6.0 


0.75 1.00 


0.88 


86.6 


9.1 


12.1 


15.2 



Table 4: Differential cross sections for the process ep — > ebbX — > ejjfiX' for the photo- 
production sample in the kinematic range Q 2 < 1 GeV 2 , 0.2 < y < 0.8, pf> 2.5 GeV, 
-0.55 < rf < l.l,pf* 1(2) > 7(6) GeV and |^' e *| < 2.5. 
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Measurement 


Experimental Errors 


Q 2 -range 
[GeV 2 ] 


Q' 2 

[GeV 2 ] 


da/dQ 2 
[pb/GeV 2 ] 


stat. 
[pb/GeV 2 ] 


syst. 
[pb/GeV 2 ] 


total 
[pb/GeV 2 ] 


2.0 5.0 
5.0 18.0 

1 Ci /"\ 1 /A/A /A 

18.0 100.0 


3.5 
9.5 
45.0 


1.55 
0.297 
0.09 1 


0.30 

0.075 

0.016 


0.22 

0.042 

0.013 


0.37 

0.086 

0.020 


rr-range 


\ogx 


da/dlogx 
[pb] 


stat. 
[pb] 


syst. 
[pb] 


total 
[pb] 


-4.5 -3.8 
-3.8 -3.1 
-3.1 -2.4 


-4.15 
-3.45 
-2.75 


6.40 
9.72 
6.68 


1.33 
l.7l 
1.44 


0.90 
1.36 
0.93 


1.60 
2.18 
l.7l 


rf -range 




da/dri^ 
[pb] 


stat. 
[pb] 


syst. 
[pb] 


total 
[pb] 


-0.75 -0.12 
-0.12 0.50 
0.50 1.15 


-0.4 
0.2 
0.8 


5.36 
7.40 
12.6 


1.42 
1.64 
1.9 


0.75 
1.03 
1.8 


1.60 
1.94 
2.6 


p^-range 
[GeV] 


ft 
[GeV] 


da/dpt 
[pb/GeV] 


stat. 
[pb/GeV] 


syst. 
[pb/GeV] 


total 
[pb/GeV] 


2.5 3.0 
3.0 3.8 
3.8 12.0 


2.8 
3.4 
6.4 


ll.3 
8.05 
0.622 


2.4 

1.39 

0.124 


1.6 

1.13 

0.087 


2.9 

1.79 

0.151 


pg-jf -range 
[GeV] 


^Breit 
rt, jet 

[GeV] 


da/dp»ff 
[pb/GeV] 


stat. 
[pb/GeV] 


syst. 
[pb/GeV] 


total 
[pb/GeV] 


6.0 8.5 
8.5 12.0 
12.0 30.0 


7.2 
10.0 
18.5 


2.20 
1.96 
0.183 


0.52 
0.31 
0.043 


0.3 1 
0.27 
0.026 


0.60 
0.42 
0.050 



Table 5: Differential cross sections for the process ep — > ebbX — > ejfxX' for the electropro- 
duction sample in the kinematic range 2 < Q 2 < 100 GeV 2 , O.l < y < 0.7, p?> 2.5 GeV, 
-0.75 < ?f < l.l5,pg e f > 6 GeV and \r] jet \ < 2.5. 
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Photoproduction 

a(ep — > efr&X — > ejj/iX') [pb] 


Electroproduction 

a(ep — > efrfrX — > ejfiX') [pb] 


Data 
FMNR 
PYTHIA 
CASCADE 


38.4 ±3.4 ±5.4 
23.81^ 
20.9 
22.6 


Data 
HVQDIS 
RAPGAP 
CASCADE 


16.3 ±2.0 ±2.3 

Q.oirJ 

6.3 
9.8 



Table 6: Measured cross sections with their statistical and systematic errors and corresponding 
predictions from NLO QCD calculations and Monte Carlo simulations in the kinematic range 
Q 2 < 1 GeV 2 , 0.2 < y < 0.8, $> 2.5 GeV, -0.55 < if < 1.1, pf* 1(2) > 7(6) GeV and 
\i] jet \ < 2.5 (photoproduction) and in the kinematic range 2 < Q 2 < 100 GeV 2 , 0.1 < y < 0.7, 
p?> 2.5 GeV, -0.75 < rf < 1.15, pg"«* > 6 GeV and |V' e '| < 2.5 (electroproduction). 
The errors for the predictions from FMNR and HVQDIS give the systematic uncertainties as 
estimated from scale variations (see text). 
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Figure 2: Distributions in photoproduction of a) the muon transverse momentum p%, b) the 
pseudo-rapidity of the muon if, c) and d) the transverse momenta p J t e 1(2) of the highest and 
the second-highest p t jets, respectively and e) the observable x° bs . Included in the figure are 
the estimated contributions of events arising from b quarks (dashed line), c quarks (dotted line) 
and light quarks (dash-dotted line). The shapes of the distributions from the different sources 
are taken from the PYTHIA Monte Carlo simulation and their relative fractions are determined 
from a fit to the two-dimensional data distribution of p T t el and 5 (see text). 
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Electroproduction 
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Figure 3: Distributions in electroproduction of a) the photon virtuality Q 2 , b) the inelasticity y, 
c) Bjorken x, d) the muon transverse momentum and e) the transverse momentum Vt^tt °f me 
selected jet in the Breit frame. Included in the figure are the estimated contributions of events 
arising from b quarks (dashed line), c quarks (dotted line) and light quarks (dash-dotted line). 
The shapes of the distributions from the different sources are taken from the RAPGAP Monte 
Carlo simulation and their relative fractions are determined from a fit to the two-dimensional 
data distribution of p r t el and 5 (see text). 
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Photoproduction 
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Figure 4: Distributions in photoproduction of a) the impact parameter 5 of the muon track and 
b) the transverse muon momentum p r t el relative to the axis of the associated jet. Included in the 
figure are the estimated contributions of events arising from b quarks (dashed line), c quarks 
(dotted line) and the light quarks (dash-dotted line). The shapes of the distributions of the 
different sources are taken from the PYTHIA Monte Carlo simulation and their relative fractions 
are determined from a fit to the two-dimensional data distribution of p T t el and 5 (see text). 
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Electroproduction 
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Figure 5: Distributions in electroproduction of a) the impact parameter 5 of the muon track and 
b) the transverse muon momentum p r t el relative to the axis of the associated jet. Included in the 
figure are the estimated contributions of events arising from b quarks (dashed line), c quarks 
(dotted line) and the light quarks (dash-dotted line). The shapes of the distributions of the dif- 
ferent sources are taken from the RAPGAP Monte Carlo simulation and their relative fractions 
are determined from a fit to the two-dimensional data distribution of p T t el and 5 (see text). 



26 



Photoproduction 
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Figure 6: Distributions in the restricted photoproduction sample of a) the impact parameter 
5 for events with p r t el > 1.2 GeV and b) the transverse muon momentum p r t el relative to the jet 
axis for tracks with impact parameter 5 > 0.01 cm. The predictions for the contributions to 
the restricted sample from b events (dashed line), c events (dotted line) and light quark events 
(dash-dotted line) , as determined from a fit to the two-dimensional distribution of p[ el and 8 in 
the full data sample (see text), are also shown. 
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Figure 7: Distributions in the restricted electroproduction sample of a) the impact parameter 
5 for events with p r t el > 1.2 GeV and b) the transverse muon momentum p r t el relative to the jet 
axis for tracks with impact parameter 5 > 0.01 cm. The predictions for the contributions to 
the restricted sample from b events (dashed line), c events (dotted line) and light quark events 
(dash-dotted line), as determined from a fit to the two-dimensional distribution of p[ el and 5 in 
the full data sample (see text), are also shown. 
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Figure 8: Differential cross sections for the photoproduction process ep — > ebbX — > ejjfiX' 
in the kinematic range Q 2 < 1 GeV 2 , 0.2 < y < 0.8, p£> 2.5 GeV, -0.55 < ^ < 1.1, 
^je 1(2) > y^-g-j Q £ y an( j |^j e t| < 2 5 cross sections are shown as functions of a) the muon 
pseudo-rapidity r] 11 , b) the muon transverse momentum pf, c) the jet transverse momentum p 3 t etl 
of the highest transverse momentum jet and d) the quantity x° bs . The inner error bars show the 
statistical error, the outer error bars represent the statistical and systematic uncertainties added 
in quadrature. The NLO QCD predictions at the parton level (dashed line) are corrected to the 
hadron level (solid line) using the PYTHIA generator. The shaded band around the hadron level 
prediction indicates the systematic uncertainties as estimated from scale variations (see text). 
Predictions from the Monte Carlo generator programs CASCADE (dotted line) and PYTHIA 
(dash-dotted line) are also shown. 
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Figure 9: Differential cross sections for the electroproduction process ep ■ 
the kinematic range 2 < Q 2 < 100 GeV 2 , 0.1 < y < 0.7, $> 2.5 GeV, 
Ptjtf > 6 GeV and \r]^ et \ < 2.5. The cross sections are shown as functions of a) the photon 
virtuality Q 2 and b) the Bjorken scaling variable x. The inner error bars show the statistical 
error, the outer error bars represent the statistical and systematic uncertainties added in quadra- 
ture. NLO QCD predictions at the parton level (dashed line) are corrected to the hadron level 
(solid line) using the RAPGAP generator. The shaded band around the hadron level prediction 
indicates the systematic uncertainty as estimated from scale variations (see text). Predictions 
from the Monte Carlo generator programs CASCADE (dotted line) and RAPGAP (dash-dotted 
line) are also shown. 



30 



Electroproduction 



73 

D 
73 



10 



ep -> ebbX -> ej^X 
• Data 

NLO QCD ® Had 

NLO QCD 

Cascade 
Rapgap 



H1 

a) 



2<Q 2 <100GeV 2 



I 



i 



-0.5 



0.5 



> 

o 
O 

B 

Q. 

'I 12 

Q. 
73 



10 



ep -> ebbX -> ejnX |_j -| 

• Data 

NLO QCD ® Had 

NLO QCD 

Cascade 

— Rapgap ^ 



i i 



2<Q <100 GeV 



10 



15 



20 



Cm 



30 



> 

0) 

O 

B 

a. 10 

Q. 
73 

C 
73 



i 



ep -» ebbX -> ej^X |_| -| 

• Data 

I NLO QCD® Had 

I NLO QCD 

Cascade 
Rapgap ' 

2<Q 2 <100 GeV 2 



I 



10 
pftGeV] 



Figure 10: Differential cross sections for the electroproduction process ep — > ebbX — > ejfiX 1 
in the kinematic range 2 < Q 2 < 100 GeV 2 , 0.1 < y < 0.7, f t > 2.5 GeV, -0.75 < if < 1.15, 
pfjef > 6 GeV and \r]^ et \ < 2.5. The cross sections are shown as functions of a) the muon 
pseudo-rapidity r/ M , b) the muon transverse momentum p% and c) the transverse momentum 
pfjef °f me leading jet in the Breit frame. The inner error bars show the statistical error, 
the outer error bars represent the statistical and systematic uncertainties added in quadrature. 
The NLO QCD predictions at the parton level (dashed line) are corrected to the hadron level 
(solid line) using the RAPGAP generator. The shaded band around the hadron level prediction 
indicates the systematic uncertainty as estimated from scale variations (see text). Predictions 
from the Monte Carlo generator programs CASCADE (dotted line) and RAPGAP (dash-dotted 
line) are also shown. 
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Figure 11: Ratio of beauty production cross section measurements at HERA to NLO QCD 
predictions. The results of this paper (solid circles and squares) are compared with ratios de- 
termined using the measurements taken from [5-7]. The photoproduction points are plotted at 
different horizontal positions for better visibility. Note that cross section definitions and kine- 
matic ranges are somewhat different for the different data points. The dotted lines indicate the 
typical theoretical error due to scale uncertainties. The theoretical prediction used to form the 
ratio with the measurement by the ZEUS Experiment [6], shown as an open square, is calcu- 
lated using the same program and parameter choices as for the prediction for this measurement 
(full square). Different parameter choices, e.g. for the modelling of the hadronisation and decay 
of the 5-hadron, lead to a variation of the prediction of ~ 10%. 
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